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Abstract 

We compute the production of quark-antiquark pairs in high energy col- 
lisions between a small and a large projectile, as in proton-nucleus colli- 
sions, in the framework of the Color Glass Condensate. We derive a gen- 
eral expression for quark pair-production, which is not fcx-factorizable. 
However, fcx-factorization is recovered in the limit of large mass pairs or 
large quark-anti-quark momenta. Our results are amenable to a simple 
interpretation and suggest how multi-parton correlations at small x can 
be quantified in high-energy proton/deuteron- nucleus collisions. 

1 Introduction 

In an accompanying paper [1], heretofore referred to as I, we developed the Color 
Glass Condensate (CGC) framework [2-4] to address the case where the parton 
densities in the projectile are small while those of the target are large. This 
may be applied for instance to the proton/deuteron- nucleus collisions that are 
currently being studied at Brookhaven's Relativistic Heavy Ion Collider (RHIC) 
and will be further studied at the Large Hadron Collider at CERN in the near 
future. In our previous paper [1], we solved the Yang-Mills equations, for two 
light-cone sources, to first order in the density of the dilute projectile and to 
all orders in the density of the target. Our computations were performed in 
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the Lorcnz/covariant gauge df^A'^ = 0. Paper I focused on developing a gen- 
eral formalism for gluon production at high energies and on phenomenological 
consequences such as the Cronin effect. In I, we derived a compact expression 
of the gauge fields produced in proton-nucleus collisions. In this paper, we will 
use this result to study the production of quark-anti-quark pairs in high energy 
proton- nucleus collisions^. 

Quark pair production in high energy hadronic collisions has prcvioTisly 
been studied in the frameworks of collinear factorization [7-9] and fc_L-factoriza- 
tion [10-12]. A comprehensive list of references for the particular applications of 
these formalisms to proton-nucleus collisions can be foimd in ref. [13]. At high 
energies, where the intrinsic transverse momenta of the partons is not negligible 
and may be on the order of the saturation scale [12], k± ^ Qs, fcx-factorization is 
a convenient formalism to discuss particle production. Phenomenological studies 
for hadron colliders have been discussed in refs. [14,15] and for heavy-ion colli- 
sions in ref. [16]. fc_L-factorization was proven for large mass pair-production at 
small X in refs. [10-12]. 

The issue of -factorization for quark production in the CGC approach 
was addressed by two of us in [17]. We showed that, in the CGC language, 
the fc^-factorization of Collins and Ellis [10] is recovered exactly^ when the 
parton densities of the projectile and the target are not too large, namely, when 
Pp/k^, Pj^/k\ ^ 1. Here Pp, p^ denote the number of parton sources per unit 
area in the projectile and target respectively and k± is a typical momentum of 
the produced parton. Our results also suggested however that fc_L -factorization 
is not robust and would fail when either or both Pp,Pj^ ^ k\. High energy 
proton-nucleus collisions are ideal to test this conjecture explicitly since the 
proton source may be taken to be dilute {pp/k^ <C 1) while the nuclear target 
is dense {pj^/k\ ^ 1). For gluon production, we (and others previously) showed 
that fc^-factorization remains valid even in proton-nucleus collisions, although 
one has to use an "unintegrated nuclear gluon distribution" which is not the 
canonical one. However, we will show in this paper that A;_L-factorization fails 
unequivocally for pair production in pA-coUisions. Nevertheless, our final result 
(eq. 48) can be expressed in terms of a small number of distribution functions 
which have a simple interpretation. 

Unlike the unintegrated distribution discussed in the case of gluon produc- 
tion, these distributions cannot all be expressed in terms of 2-point correla- 
tors of adjoint Wilson lines. In addition to a term containing the correlator 
{U{xi_)U\y^}), one has a term proportional to {U{xi_)t''U\y ^^f^U'^^ix'^) , 
i.e. a product of two Wilson lines in the fundamental representation (denoted 
U) and an adjoint Wilson line, and a third term proportional to the product of 
four fundamental Wilson lines: {U{xi_)t''U'^{y^^UHy'^)t''U'^{x'J). The trans- 
verse coordinates here represent those of the quark and the anti-quark in the 
pair production amplitude and in the complex conjugate amplitude. For very 

^Note that this problem has also been investigated numerically by Dietrich using a toy 
model for the background classical color field [5,6]. 

^The results of Collins and Ellis are equivalent to those of Catani, Ciafaloni and Haut- 
mann [11]. 
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large mass pairs and for large momenta of the pair, the transverse separation 
of the pair can be neglected, and one can show, through an identity, that the 
A;_L-factorized result in terms of the two point function of adjoint Wilson lines 
is recovered. 

All the above Wilson lines are computed in the presence of the classical color 
field of the nucleus and therefore depend on the source p^. The brackets (• • • ) 

here denote the weighted quantum mechanical average of an operator in the 
background field of the proton and the nucleus, which can be expressed as 

(O) = j [Dpp][DpM^o,PpWAx^,Pa] o[pp,pA ■ (1) 

This averaging procedure is essential in order to restore gauge invariance since 
0[pp, p^] is computed in a particular gauge. More importantly, it is also through 
Wp and Wj^ that quantum effects, due to evolution of the light cone wave func- 
tions of the target and projectile with x, are incorporated in the study of pair 
production at high energies. The arguments and .Tq denote the scale in x 
separating the large- a; static sources from the small- a; dynamical fields. In the 
McLerran-Venugopalan model, the functional Wj^ that describes the distribu- 
tion of color sources in the nucleus is a Gaussian in [18 21]. In general, 
this Gaussian is best interpreted as the initial condition for a non-trivial evo- 
lution of W^^[a;o with Xq. This evolution is described by a Wilson renor- 
malization group equation (often called the JIMWLK equation [22-33]). This 
evolution equation has been solved recently numerically by Rummukainen and 
Weigert [34] . One is therefore in a position to make predictions for high energy 
pair-production in hadronic collisions that arc sensitive to not only the "unin- 
tegrated gluon distribution" but to more subtle correlations among partons in 
hadrons at high energies. Conversely, comparisons with pair production exper- 
iments at high energies can test, with a high degree of precision, whether the 
JIMWLK equations, or equivalently the Color Glass Condensate, is the right 
effective theory of high energy QCD. 

This paper is organized as follows. In section 2, wc briefly review the solution 
to the Yang-Mills equations derived in I and re-express it in a way convenient for 
our discussion. We next proceed to the computation of the quark pair produc- 
tion amplitude in section 3. This amplitude consists of two sets of terms which 
we name "regular" and "singular" terms. The former correspond to the follow- 
ing cases: a) the gluon produces the pair before interacting with the nucleus - 
the pair therefore scatters on the nucleus on its way out, and b) the gluon inter- 
acts with the nucleus and produces the pair after the encounter. The singular 
case corresponds to the situation where the pair is produced inside the nucleus 
and re- interacts on the way out. Naively, one would not expect this term to 
contribute at very high energies, due to the brevity of the encounter. Neverthe- 
less, in the Lorenz gauge that we choose to work in, such a term does exist and 
is important in obtaining our final result. Adding both sets of terms together, 
we obtain the time-ordered amplitude for pair production. For completeness, 
and to confirm our result for the time ordered amplitude, we also compute the 
retarded pair production amplitude, which, to this order, carries an identical 
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physical content albeit a different interpretation. The computation is presented 
in appendix B. In section 4, we compute the pair production cross-section. Our 
final result (eq. 48) is not fc_L-factorizable in the same manner as the result for 
the gluon production cross-section is. The various terms in our result however 
have a simple interpretation and for large mass or large momentum pairs, one 
does recover the A;_L-factorized result - as expected from our derivation in [17]. 
We finally address the case of the single quark production cross-section in sec- 
tion 5. Single quark distributions were recently studied by Tuchin [35] in the 
context of a Glauber model of independent scatterings. In contrast to our mo- 
mentum space approach, Tuchin's results are formulated entirely in coordinate 
space and therefore do not address the issue of A;_L-factorization. Further, our 
results are more general, being valid for non-Gaussian correlations as well. For 
the case of a Gaussian distribution of color sources, we provide in appendix A 
explicit expressions for all the correlators that appear in our results for both 
pair production and inclusive single quark production. We end with a summary 
of our results and an outlook on future work on these questions. 



2 The Yang-Mills gauge field produced in pA 
collisions 

Our first step in calculating the pair production amplitude is to find solutions of 
the classical Yang-Mills equations in the presence of proton and nuclear sources 
of color charge. Since we want to compute quark production at the lowest order 
in the source pp that describes the proton and to all orders in the source 
that describes the nucleus, we need to compute the gauge field up to the same 
order. The Yang-Mills equations, to this required order, were solved in covariant 
gauge in I. We will here just quote the results and re-express them in a manner 
convenient for the computation of the pair production amplitude in the next 
section. Following I, we denote the field of the nucleus alone, the field of 
the proton alone and (without any subscript) the total gauge field at order 
1 in the proton source and to all orders in the nuclear source. 

The gauge field created by the nucleus alone is given in coordinate space by 
the expression 

A>^^{x) = -95'^-5{x+)^p,{x^). (2) 
Its Fourier transform, which will appear repeatedly as we proceed, is 

A^M = 2n95^-5{q+)P^^. (3) 

Note that this field is linear in the source p^: it is indeed a well known fact 

that, in covariant gauge, the classical color field created by a single projectile, 
no matter how dense, is linear in the classical color source. We denote as 
the first correction to this field due to the proton source, which is of order one 
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in Pp. This correction was derived in the covariant gauge in I (see section 3.5 of 
I)- 

Our expression for Ai^ has both singular terms (proportional to 5{x'^)) and 
regular terms. It will be convenient for our later discussion to group them 
accordingly, and express as 



A^{q)=A>i,^{q) + 5^-A^^^{q). 
Our expression for A(4g is 



(4) 



^9 



C^iqMi) [U{k2±) - (27r)25(fe2±)] 



q^ + iq+ej {2tt)^ 

+C^,e,{l) [V{k2^) - {2nrS{k2^)] ]^^^ . (5) 



+ 



In this formula, the first term is the color field of the proton alone, given by 
eq. (3) with replaced by pp. In the second term, feij_ is the momentum 
coming from the proton and k2±, defined as k2± = q^ — ki±, is the momentum 
coming from the nucleus. The 4-vectors C^{q, ki±) and C^^^egl^) '"^^^ given by 
the following relations: 



C+{q,k,^)^ 



-^ii- , C-{q,k,^) ^ , Cl{q,k^^) ^ -2k\ 

q +ie " q+ ^ 



C+res{q) = 2q+ , C- ,,g(9) = 2^- - ^ , C;,,,g(g) = 2gV (6) 

U and V are the Fourier transforms of Wilson lines in the adjoint representation 
of SUiN): 



U{k^) 
with 



d^x^e'-'^^ '^^Uix^) , V{k^)= d^x_Le"'^ '^^V{x_L) 



(7) 



U{x±) = P+ exp 
V{x±) = V+ exp 



^9 



/-hoo 
dz+A-{z+,x^)-T 
-OO 

/+00 
dz+A-{z+,x^)-T 
-OO 



(8) 



where the T'^ are the generators of the adjoint representation of SU {N) and P+ 
denotes a "time ordering" along the z~^ axis. The subscript 'reg' on indicates 

that the corresponding term of A^' docs not contain any 5{x^) when expressed 
in coordinate space. The relationship^ between reg (9) the vector Ci^{q) 

*The terms responsible for a i5(a;+) in Ai^(x) are those for which C^^jq^ does not go 
to zero when q~ — > oo. It is immediate to verify that this happens only for C~. C^^reg is 
obtained from C~ the piece responsible for this non-zero limit. 
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introduced in the section 3.5 of I is^: 



C^^;,reg(9)=^?(^(9) + 5''-^. (9) 



The "singular" term reads: 

J2 



47ng(«) - / ^ [nk..) - {2.f5{k,^)] . (10) 



It will in fact be more convenient later to have this part of the field in coordinate 
space: 

where V{x^ ,—oo; x±) denotes an incomplete Wilson hue: 



A;,^^{x)=i^iA-{x)-T)V{x+-oo;x^)9{x-)-^Pp{x^) , (11) 



V{x~^, — oo; x±) = V+ exp 



r dz+A-{z+,x^)-T 

J — OO 



(12) 



In the case of gluon production I, we have seen that the Wilson line V drops out 
of the gluon production amplitude. This result was to be expected from the fact 
that there arc gauges in which such a Wilson line, characterized by an unusual 
factor 1/2 in the exponent, does not appear at all. Naturally, we expect the 
same for quark-antiquark production, and one of our tasks in the next section 
will be to exhibit the mechanism by which this Wilson line drops out of the pair 
production amplitude. 

3 Pair production amplitude 

3.1 Generalities 

Once we know the time-dependent classical gauge field created in the collision 
of the two projectiles, it is straightforward to calculate the production of quark- 
antiquark pairs in the collision. One can simply forgot about the proj(X'tiles 
themselves, and consider only the classical field. The probability to produce 
exactly one qq pair^ in the collision can be expressed in terms of the time- 
ordered quark propagator in the presence of the classical field^ as 



^Cjj and are related to the well known Lipatov vertex [36,37] via the relation = 

+ \ . The properties of these arc explored further in paper I. 

^If one uses instead the retarded quark-propagator, one obtains the average number Nqq 
of produced pairs in the collision. At leading order in pp, at most one pair can be produced 
per collision, so that one expects Pi = Nqq. We check this property explicitly in appendix B. 

^This formula for Pi is correct at leading order in pp, but is incomplete in general. Indeed, 
it should also contain the prefactor |{0in|0out)|^, which is mandatory for unitaxity to be 
preserved (see [38] for more details). This prefactor only affects higher order corrections in 
Pv 
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where the argument [pp,/3^] indicates that this is the production probabihty in 
one particular configuration of the color sources. In order to turn this probability 
into a cross-section, one must average over the initial classical sources pp and 
with the weight Wp[pp]W^[/9^] and one must integrate over all the impact 
parameters 6, to obtain, 

ai =1 dPb J [Dpp][DpjWMWMP,[pp,pJ . (14) 

In eq. (13), the time-ordered amplitude is related to the quark propagator 
as follows: 

M^{q, p) = u{q)T^ {q, -p)v{p) , (15) 

where (g, —p) is the interacting part of the Feynman quark propagator in the 
presence of the classical field, with its external legs amputated. The arguments 
q and —p are respectively the outgoing and the incoming 4-momenta on this 
quark propagator. 

3.2 Regular terms 

Since we want the pair production amplitude to first order in the source pp and 
to all orders in , we must insert the total field exactly once on the quark 
propagator, and we can insert the field A'^ an arbitrary number of times since 
it does not change the power counting in pp. At this point, it is convenient to 
study separately the terms in A^ which do not contain a 6{x~^) from those that 
contain a 5(x+). 

Let us begin with the former. The terms that do not contain a S{x'^) are the 
ones in eq. (5) that we called A{4g- Physically, these are terms for which the pair 
is produced outside the nucleus; more precisely, the point where A^^ is inserted 
on the quark line must be outside the nucleus. Regular terms are represented 
by the four diagrams of figure 1. In order to calculate these diagrams, we need 
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(a) (b) 




Figure 1: Regular terms in the time-ordered pair production amplitude. The 
gluon line terminated by a cross denotes a classical field insertion. The black 
dot denotes multiple insertions of the field A'^ (with at least one insertion). 
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the expression for the multiple insertions of the field on the time-ordered 
quark propagator, denoted by a black dot in the diagrams. This reads [38,39]: 



p+k 
m > — 



27r(5(fc+)7+sign(p+) / d?x^e'''^ ''^ \w'^''^p^\xi_) - I 



(16) 



where [/ is a Wilson line in the fundamental representation of SU{N), defined 
by 



U{x±) = V+ exp 



dz~^A^ (2+,x_l) • t 



(17) 



The t"' are the generators of the fundamental representation of SU{N). The 
contributions of the diagrams of figure 1 to the time-ordered pair production 
amplitude are obtained using the free Feynman (time-ordered) quark propagator 
in the intermediate quark lines: 



1^ + m 



' — rn? + ie 

The contribution of the diagram (a) to Al^ is simply given by 



(18) 



(19) 



Since the first term of eq. (5) brings a S{p~ + q~) with it, this term cannot 
contribute to the production of a pair on shell, for which we must have p~ > 
and g~ > 0. (This first term would correspond to the production of a pair from 
the proton alone, which is kinematically forbidden). Using also u{q)(;^+4)v{p) = 
0, we can write the contribution from diagram (a) more explicitly as 



M^^\q,p)=g^ j 



11. Pp,a (fci±) 



u{q)^u{p + q,ki±)t^v{p) 



u{q)^^t^v{p) 



p+ + q 



?^[V(»^x)-l|4. 

(20) 



The contribution of the diagram (b) reads: 



M 



q,p) = -J ^27rS{k+)[-tgA^,lip + q-k)] J d'x^e 
xu{qhy.eSl{-p + k)j+ \W{x^) - l] v{p) 



ik± -xj 



(21) 



where we have already used the fact that the quark propagator on which the 
eikonal interaction with the nucleus is inserted carries a negative energy —p'^ < 
0. The integral is trivial because of the 6{k^). The k~ integral can be 
performed by using the theorem of residues. The propagator 5° (— p + k) has a 
pole above the real axis in the k~ plane. So do the denominators {p + q — fc)^ 
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(the prefactor in the second term of eq. (5)) and p~ + q~ — k~ + ie (in C^). 
Therefore, we see that the second hne in the field cannot contribute because 
it has all its poles on the same side of the real axis . The only term from Ai^^^ 
that contributes is the first one. namely, the field of the proton alone. Finally, 
using the variable feij, = p_L + <7j^ — fej. instead of k^_, we have 



uiqh-f'ii -h+ m)-i+[U\xi) - l]v{p) 



2p+q- + (q_L - ki±f + 



(22) 



Similarly for the diagram (c), the poles in the variable k are such that only 



the first term in Af^^^ contributes. We obtain 



u{q)^+[U{xi_) - l](-y+|^i +m)7-ri;(p) 



(23) 



2p- q+ + {p ^ - kn_Y + m? 
By the same method, we can write the contribution of the diagram (d) as follows: 

M^f{q,p)=g'l "12^)2 'I2)'2 '"i^l"-^ jd-x^d^y^e^'^.-^.e^iv.^^^M-v. 

^ Xlh^[U{x±)-^\ii -t + m-h-f'ii m)j+[U^y^)-l]v{p) 
2?'+[(9_L - fc-i)^ + m^] + 2«+[(q_L - fe-L - fei±)^ + ^2] 

(24) 

In this expression, x± and yj^ are the transverse coordinates of the quark and the 
antiquark respectively, kij_ is the momentum flowing from the proton, k±_ the 
momentum flowing from the nucleus to the quark line, and pj^ + qj^ — k± — ki± 
the momentum flowing from the nucleus on the antiquark line. 

One can sec at this point that the sum of the regular terms'' contains the 
Wilson line V. As we shall see shortly, it is cancelled by the contribution coming 
from the singular terms. 



3.3 Singuleir terms 

We must now evaluate the contribution to the pair production amplitude coming 
from the held A'^^^^. Since this field is proportional to d{x^) in coordinate 
space, its contribution is a term where the pair is produced inside the nucleus. 
Technically, this means that in order to correctly compute this term, we must 

*For this conclusion to be valid, it was crucial to have isolated the singular terms in the 
gauge field from the regular ones. Indeed, the singular terms have a large k~ behavior which 
prevent one from closing the integration axis by the contour at infinity. 

^In fa<;t, only the diagram (a) contains the Wilson line V. For the diagrams (b), (c) and 
(d), only the term of the total field A'' contributes, for kinematical reasons. 



9 



go to coordinate space and regularize the 5{x~^) by giving a small width to the 
nucleus, 

5{x+) ^ 5,{x+) , (25) 

where ^e(x+) is a regular function (whose support is [0, e]), which becomes a 

6{x^) when g goes to zero. We need not specify this function further because the 
results do not depend on the precise choice of the regularization. The singular 
contribution is depicted in figure 2. The field A'^-^^ is inserted on the quark line 
at the 'time' the quark and the anti-quark then rescatter on the field 
of the nucleus between x~^ and e. This contribution to the amplitude can be 













A . Womoo 
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e 


"sing ^ooQooo 




H 

Aa 





Figure 2: Singular term in the time-ordered pair production amplitude. The 
gluon line terminated by a cross denotes a classical field insertion. 



written in coordinate space as, 

A^';"^(g,P)= J d^xe'^P+''>''u{q)Ui+oo,x+;x^) 

where we have introduced the incomplete Wilson line 



U{+oo, x^;x±) = ■P-i- exp 



r+<x) 

/ dz+A-{z+,x±)-t 

Jx+ 



(26) 



(27) 



Note that here the fundamental Wilson line and its complex conjugate are eval- 
uated at the same point x_\_. This is because in the limit e ^ 0, at the end 
of the interaction, the components of the pair will not have had sufficient time 
to separate from each other. Note further that the four possible contributions 
(analogue to the four regular terms) are already summed here by replacing U —1 
by U and C/t - 1 by . 

In order to simplify this expression, we need the following algebraic identity: 

U{+oo, x+; XA_)t''U\+oo, x+; x_l) = t*?7''"(+oo, x+; x_l) . (28) 

Using also the explicit form eq. (11) of the singular field, we can write the 
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singular contribution to the amplitude as^°, 

x|i| j dx+[U{+<x^,x+;xA_){A-{x)-T)V{x+,-cx;xA_)\^}^ 
xu{q)^+t''e{x-)^Pp,b{x±)v{p) . (29) 

Finally we need the formula 

+00 

i| j dx+U{+oo,x+;x_L){A-{x)-T)V{x+,-oo;x±) = U{x_L)-V{x_L) , (30) 



which we encountered previously in I. The remaining Fourier integral in eq. (29) 
can be performed exphcitly, and wc obtain 

x^^^^?^$^[''(xj-a(xj]... (31) 

We see that, as anticipated, the term in V in A^^"^s cancels the term in V we 
found previously in 

3.4 Complete time-ordered amplitude 

We must now combine the terms A^(^">'''^''') and A^""s order to do so, we 
shall use the following identities: 

l~v{p) = ^'y~^- m)'y+v{p) , (32) 

which are simple consequences of the Dirac equation obeyed by the free spinors. 
The second trick is to introduce dummy variables fcj, and yj^ via the identity 



^"We have used the fact that the support of the a;+ integral is < a;+ < e in order to 
approximate exp(i(p~ + q~)x^) 1. 
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so that all the terms have the same integration variables. Putting all these terms 
together, we obtain the following expression for the complete amplitude^^: 



gifej_-a!± ^i{p±+q±-k±-ki±)-yj_ 



+u{q) 



2p+[{q^ - fei)2 + m2] + 2q+[iq^ - k± - feix)^ + m^] 



t''v{p)U''%x^)} . (34) 



{p + q? P+ + q+_ 

This is our final formula for the time ordered pair production amplitude. 

3.5 Physical interpretation and limits of eq. (34) 

The interpretation of the two terms in this amplitude can be made more trans- 
parent if we note that 



u{q) 



(p^{p + q,ki±) 7+ 



(p + q)"^ p^ + q 



+ 



fviP) = ^(g) (^^g)2 ^ Mp) , (35) 



where C^(p+g, ki±) is the efi'ective Lipatov vertex for the production of a gluon 
of momentum p + q via the fusion of two gluons of momenta fci and p + q — ki. 

We can then interpret the first term of eq. (34) as the production of a qq pair 
by a gluon from the proton before it goes through the nucleus. It is therefore 
the qq state that interacts with the nucleus, and hence the two Wilson lines 
in the fundamental representation at two different transverse coordinates. The 
second term in eq. (34) corresponds to a gluon from the proton going through 
the nucleus before producing the qq pair. We recognize the Lipatov vertex that 
is known to appear in the gluon produc;tion amplitude in pA collisions (see the 
discussion in I) and the adjoint Wilson line U describes the rescatterings of 
this gluon while it goes through the nucleus before it finally produces the pair 
outside the nucleus. 

It is straightforward to verify that eq. (34) agrees with the expression that 
we obtained in [17] in the limit of low nuclear density (see eqs. (15), (19) and 
(20) of [17]), namely, when we expand all the Wilson lines of eq. (34) to first 
order in the source . 

The amplitude in eq. (34) satisfies another important property: there is a 
cancellation between the two terms of eq. (34) in the limit ki± — > 0. Indeed, for 
ki± = 0, the curly bracket of eq. (34) becomes 

r . _ u{q)^+[U{x^)t-m{yMp) _ uiqh+[t''U"'{x±Mp) 
^'"i- p++q+ p++q+ ■ ^"^^^ 



^^In order to obtain this expression, we have also used a cancellation between the —f^ terms 
in Uf'U^ - t" and in t^[U - l]ba- 
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Since this expression does not depend on fc_L any more, the integration over fe^ 
in eq. (34) produces a 6{x± —y±)- We can therefore replace U{x±)t°'U^yj_) in 
the first term by t^U^"'{x±) (thanks to eq. (28)) - which proves the announced 
cancellation. This property ensures collinear factorizability on the proton side^^ , 
by making the ki± integral at most logarithmically singular. 



4 Pair production cross-section 

4.1 Average over the color sources 

Before we square the amplitude in order to obtain the pair production cross- 
section, it is convenient to write it in a more compact form by introducing some 
shorthand notation. We shall write the amplitude as 



^ \ 2 /(^^fel-L C^^fe-L Pp,a(fel_L) f,2 ,2 



ik± x± ^i{p±+q±-k±-ki±)-yj^ 



xuiq)^T,gik,^,k^)[Uix^)t^U\y^)]+Tg{ki^)[t''U'''{x^)] jv{p) . 
where we denote^^ 



(37) 



Tgg{kl±,k^) = 



2p+[{q^~k^)^ + to2] + 2g+[(q^-fei-feii)2 + m2] 



Squaring the amplitude, and averaging over the configurations of the color 
s()ur( (>ri f>i, and /) ^. \\(^ luwc llie following expression for the differential proba- 

^^This property has been assumed and exploited in several recent studies of pA collisions [40- 
44]. 

^'^Thc momenta p and q of the produced particles have not been listed among the arguments 
of these objects in order to make the equations more compact. 
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bility to produce a pair in a collision at impact parameter 6: 



1.2 y2 



gi(fci-xj.-fcl-a!'j_)gi(p_L+qj_-fcj_-fcij.)-y_Lg-i(pj_+ijj_-fc^-fc'u_)-y^ 



x±..x',y^.y. 



x{trd[(^+m)Tqq-(|/-m)7%V]trc(L^(a;±)i"C/t(y^)C/(yl)r'[/t(a;l) 
+trd [(^+m)T,,-(2^-m)7°Tty]tre (c/(a;^)rC/t(y^)t6'c/ta';.'(^'^) 

+trd [(^+m)Tg(2/-m)7%V] t^c lt'>U'''{x^)U{v'^)e' U\x'^) 

+trd[«+m)T<,(y-m)7"rJV]tre (t^U'"'{x^)t''' {x'S) } 



(39) 



In this formula, trd denotes a trace over the Dirac indices (we sum over all the 
spin states of the produced fermions), trc is a trace over the color indices (we 
sum over the color of the produced quark and antiquark). We have omitted the 
arguments in Tqq and Tg, and T'^^, Tg denote the same quantities with kij_, k±_ 
replaced by k'ij_,k'j_. We can already sec that, in general, this probability 
depends on 2-, 3- and 4-point correlators of Wilson lines. As in I, we can 
introduce the proton unintegrated gluon distribution as follows: 



5^(pp,a(A;i±)p^,„,(fe'l_L)) = 



Trd. 



ki± + fe'u 



■ (40) 



The integration over X ±_ runs over the transverse profile of the proton, = 
N"^ — 1 is the dimension of the adjoint representation of SU (N), and dipp/d!^X± 
is the number of ghions per unit area and unit of transverse momentum in the 
proton (namely, the proton non- integrated gluon distribution per unit area). As 
long as we are interested in transverse momenta which are large compared to 
the inverse proton size, we can approximate feij, « k'n^ and write 

9^ {ppAk^MA^'^^)) = '^/..J e^C'^^-'.^)-^^^!^ . (41) 



We will also neglect the difference ki± — k[j_ in the other factors of eq. (39). 

Let us now turn to the various correlators of Wilson lines that appear in 
eq. (39). In eq. (39), the coordinates x±,x'j_,y^,y'j^ are relative to the center 
of the proton. Shifting all of them by the impact parameter b makes them 
relative to the center of the nucleus, at the expense of introducing an extra 
factor exp(i(fei_L — k[j_) ■ b) in eq. (39). A generic property of these correlators 
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of Wilson lines is that, for a large nucleus, they arc approximately invariant 
by translation in the transverse plane: the violations of translation invariance 
arise when the separation between two coordinates become comparable to the 

nuclear radius. 

The simplest term in eq. (39) is the fourth one, involving a correlator between 
two Wilson lines in the adjoint representation. By analogy with eq. (41), we 
can define an unintegrated gluon distribution dip^^^ / (P X ±_ for the nucleus as 
follows: 

§aa' j e*('^-^-^-''-^-^'^Hrc(t''f/'"'(a;_L)t'''f/1'"''''(a;'_L) 

Xj_ ,£C^ 



Again, we have neglected l± — Z'l in these manipulations, because this difference 
is of the order of the inverse of the radius of the nucleus or smaller. This equation 
can be seen as a definition of the unintegrated gluon distribution in the nucleus, 
that agrees at leading order in the source with what we have used for the 
proton. It is a natural object to use in this problem because it absorbs all 
the rescattering effects of the gluon on the nucleus. In fact, the same quantity 
already appeared for gluon production in pA collisions"'^''. Note that this is not 
the expectation value of the canonical number operator (except at leading order 
in the source p^). Up to some trivial factors, this object is the square of the 
scattering amplitude of a gluon on the nucleus. 

In order to deal with the 3-point correlator, it is convenient to define a 
function dipf'9/d'^X± such that: 

S' 



(43) 



27rll J ^ d^X± 



e 

Note that by construction, we have 

dipf'S{k^,l±-kr,l±\X^) d<p<^'a{l^\X^) 



I 



d?Xi_ d'^Xi 



(44) 



Indeed, integrating over fe^ forces the coordinates x±_ and y ^ of the components 
of the qq pair to become equal. When this occurs, the qq pair in an octet state 

^*<^S'S is identical to the (fj^ introduced in section 4.1 of I. The superscript g,g has been 
introduced for reasons that will become obvious later. 
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is effectively indistinguishable from a gluon. Mathematically, this statement is 
a consequence of the identity in eq. (28). One can similarly write the 4-point 
correlator as follows 



6"- I e 



i(k±-x±-k'^-x'j^)^i{l±-k±)-yj^^-i{l'^-k'j^)-y'^ 



xtre (C/(xi)i«C/t(y^)f/(yl)t«'[7t(a;'^) 



(45) 



Note further that this definition is such that 

d^f'i^{k±,l ± -k±;k'^,l±- k'JX^) 



I 



I 



(PX^ 

d^f'3{k^,l^ - k^;l^\X^) #r(«i|Xi) 



d^X_L 



d^Xj 



(46) 



In other words, collapsing one qq pair to a single point restores the 3-point 

correlator, and collapsing each pair to a point restores the 2-point correlator. 

In terms of the functions we have just introduced, we can write the differ- 
ential production probability as follows: 



dPiib) 



d^p^d^Qj^dypdyq Sir^d^ 



J S{pj_ + qj_-ki±-k2±) 



ki±,k[^,k2±_ 



I 



^^ik,^-k[^)iX^-Y^+b) d>fip{ki±\X^) 



X,,Ya 



X < / trd 

^k^,k'^ 



(^+m)Tg,-(^/-m)7°rJ,V 



d?X^ 

rfyf -^^'(fc^, fc2±-fc±; fcl, fc2±-fcY|y±) 



trd 



fe', 



+trd 



(^+m)T,(2/-m)7°Tt^7°" 

(^+m)Tg(2/-m)7°Tt'70 



#f'g(fcY,fc2±-fcl;fc2±|y±) 

dy>r(fc2±|y±) ] 

d'x^ { • 



(47) 



At this point, it is straightforward to integrate over the impact parameter in or- 
der to turn the probability Pi (6) into a cross-section. This integration produces 
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a 5{ki± — k'ljj, which makes the integration over k'ij_ triviaL Then, integrat- 
ing over the coordinates X± and Y± can be done independently in the proton 
and in the nucleus, which restores the <fp and the t/?^ functions for the whole 
projectiles rather than per unit area. Our final expression for the differential 
cross-section of quark pairs produced in proton-nucleus collisions is 



da 



atN 



<Ppj_d'^qj_dypdyq Sn'^d^ 



Hp± +Q±- ki± - k2±) 



"'1±'*2_L 




{i+m)Tg^(^-m)^Xh''l vfHku k2±-kr, k'^, k2±-k'^) 

m)T,,-(y-m)7"Tj7"l ^f'0{k±, k2±-k±; k2±) 



+m)Tg(2/-m)70Tj,-70 ipf''{k^, fe2±-fe±; fe2±) 



+trd 



(48) 



This expression for pair production is the main result of this paper. Note that 
since feij_ = k[j_, the abbreviations Tgq, T^^, Tg are related to eqs. (38) as follows: 



<-qq 



Tqq{kl^,kl_) 



T'- 

QQ 



T<,,-(fei±,fel) 



Tg = T,(fel^) . 



(49) 



The Dirac traces in eq. 48 can be computed quite straightforwardly using a 
symbolic manipulation program like FORM. We will give explicit expressions 
for these in a future paper where we explore the phenomenological consequences 
of our result. 



4.2 Leading twist approximation 

The unintegrated distributions that we encountered in the previous section con- 
tain contributions from all twists and are in general quite complicated. As 
discussed briefly in the introduction, they can be obtained from numerical so- 
lutions of the JIMWLK equations [34]. For Gaussian correlations, these can 
be computed analytically, and an explicit derivation is given in appendix A. 
Here, we will discuss only the leading twist approximation, which has a simple 
physical interpretation and is useful in order to compare eq. (48) with the result 
obtained for pp collisions. 

At leading order in the density /i^ that describes the nucleus, it is easy to 
derive a closed form expression for the functions ip^^^, f^/^^ and {pi^''i'i intro- 
duced in the previous section. The only ingredient needed for this leading order 
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calculation is the average of two 's, namely, 



(50) 



Prom the definitions given in the previous section, one then obtains the 
following relations: 



¥'r'(fe±,fe2±-fe±|A;2±) 



„2 



»'2_L 



(27r)2 [^(fei) + 5(fc±-fe2i)] 



(^f •««(fei, fe2± - fe±|feY, fe2± - feY) « 7rJ?2 



C 



+^ {S{k± - fe2±)5(fel) + 6{k^)d{k'^ - fe2^)) 



(51) 



In the derivation of these relations, we have assumed that the nuclei have a large 
transverse area ttB? and that their density is nearly uniform in the transverse 
plane. If one substitutes these leading twist approximations in eq. (48), it is 
straightforward to re-obtain the expression of the cross-section we derived in a 
previous paper for collisions (see eqs. (42), (B5), (B6) and (B7) of ref. [17]). 

Note also that the leading twist expressions for (/^^'^^ and ip'M't^ have a simple 
interpretation. At this order, the probe (gluon or qq pair) interacts with the 
nucleus by a single gluon exchange. The two delta functions in </5^''' correspond 
to the gluon being attached to the antiquark line ((5(fc±) - there is no momentum 
flow from the nucleus to the quark line) or to the quark line {5{k2i_ — fc±) ^ 
all the momentum from the nucleus flows on the quark line). An identical 
interpretation holds for the four terms of 



4.3 A;j_-factorization for pair production ? 

From cq. (48), one sees immediately that in general we cannot write the pair 
production cross-section for pA collisions in a fcj^-factorized manner, at least 
not in the usual sense. The closest to fcj_-factorization we can achieve for pA 
collisions is precisely cq. (48), where we sec that we need three different func- 
tions in order to describe the nucleus. The physical reason for this is that a 
quark-antiquark pair is a composite object: the amplitude for producing a pair 
depends on how much momentum the nucleus gives to the quark and how much 
momentum it gives to the antiquark (in case the pair is formed before going 
through the nucleus), and not just on the total momentum taken to the nu- 
cleus. This also explains why a true fc^ -factorization is possible when we treat 
the problem of pair production at leading order in the source . Indeed, at this 
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order, the pair can only exchange one gluon with the nucleus, which is either 
attached to the quark or to the antiquark, and a single momentum variable 
completely characterizes this exchange. 

From this discussion, one can also foresee in what limit a true k± factoriza- 
tion would be possible. From the integral relations of eqs. (44) and (46), we 
indeed see that if we can neglect the dependence on fej_ in Tqq{ki±, k±), then 
the integrations over k± and k'^ act only on the functions ipi^^i^ and if'^^'^ and 
simply give ^^'^ which we can now factor out. Such an approximation is valid 
when the momenta p_L,q_L of the produced quark and antiquark are large, or 
when the mass m, of the produced quark is large. In this case, one (effectively 
produces a pair which has a small transverse size and behaves like a single gluon. 



5 Single quark cross-section 

We shall now compute the single quark inclusive cross-section which is phe- 
nomenologically interesting in its own right. It is obtained from the pair produc- 
tion production cross-section given in eq. (48) by integrating over the transverse 
momentum and rapidity Dp of the antiquark. Integrating over _pj_ is trivial 
thanks to the delta function 5{p_i_ + q± — fcix — fe2±), and simply amounts to 
replacing ki± + k2± — q± in the rest of the expression. 

There is an interesting simplification of the first term, involving the 4-point 
function ipi^'^^. One notices that Tqq (see the first of eqs. (38)) does not depend 
on p^. Moreover, since Tqq has a 7"*" on the right, the product Tqq(]^—m)'^'^Tqqj° 
does not depend on either (because 7+ (5/— m)7+ = 2p+7+ does not depend 
on pj^). Therefore, this object does not depend on k2± after the substitution 
Pj^ — > fei^ -I- k2± — q±- As a consequence, the integration over k2± can be 
performed directly for this term, since we have 



^(27r)25(fe^-fel) J e''-<----'^hv,{u{x^)U\x'J) . {^2) 



Thus we see that the 4-point function reduces to a much simpler 2-point func- 
tion. It is then natural to define, by analogy with eq. (42), the following function: 

/ e^'=-(----'-)tr, (p{x^)U\x'^)) ^ -^^if{k^) . (53) 

The normalization has been chosen such that this definition becomes equivalent 
to eq. (41) in the limit where the source which describes the nucleus is a weak 
source. Because Tg depends explicitly on pj^, it becomes A;2_L-dependent after 
the substitution pj^ — > fcij, + k2± — q±- For this reason, no such simplification 
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occurs for the last three terms of eq. (48). We can now write the single quark 
production cross-section as follows: 



da„ 



aiN f dp+ 




(^+m)T,g-(feii,fc2±)(2/-m)7"rL(feii,fe2±)7°J^'Pr(fe2±) 
\4+m)Tgg{ki^, fe±)(2^-m)70Tt(fei^)70] ipf'^{k^, k2±-k±; k2±) 

{4+m)Tg{ki^)(^-m)j°T^g{ki^, k^)j°] vf\k^, k2±-kr, k2±) 



+trd 



{4+m)Tg{ki^)(^-m)j°T^„{kii_)j° 



<f'^'{k2A_) 



■ipp{k 



(54) 



In the first term, we have re-labeled k± as k2± in order to be able to factor out 
the term 1/^2.1 ■ Note also that pj^ should be understood as feij, + k2± — q± 
in this formula. Obviously, we cannot factor out in eq. (54) a common function 
(p^ that would describe the gluon content of the nucleus. This of course means 
that the single quark production cross-section, even if simpler than the pair 
production cross-section, is still not fc_L-factorizable. 

Recently, in ref. [35], Tuchin considered the inclusive single quark distribu- 
tion in proton- nucleus collisions for a model where the collisions on the successive 
nucleons are independent. This is equivalent to the case of Gaussian correla- 
tions in the MV-model. Tuchin's discussion is in the light-cone gauge = 0, 
and is expressed as the convolution of a) the square of the amplitude for the 
emission of a gluon g by a valence quark of the proton b) the square of the 
amplitude for the splitting of this gluon into a qq pair and c) the rescatterings 
of the qvgqq system^^ on the nucleus. Since his results are expressed entirely 
in coordinate space. Tuchin does not address the issue of fc^-factorization. Fur- 
ther, since our results, in contrast, are entirely formulated in momentum space, 
a detailed comparison of his result with the Gaussian correlation limit of our 
general expression is not at present feasible. Neverthek^s. several features of 
his result can be identified with terms appearing in our result. 



6 Conclusions 

We have computed in this paper the production of quark pairs in high energy 
proton-nucleus collisions in the Color Glass Condensate approach. We find that 

^^In our expression, the valence quark of the proton that emitted the gluon does not appear 
in the correlators that describe the rescatterings on the nucleus. In our language, the rescat- 
terings of this valence quark correspond to gauge rotations of the color current of the proton, 
and are taJjen care of by the Lipatov vertex. 
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the pair production cross-section, unlike the inchisive gluon production cross- 
section discussed in I, is not fc_L-factorizable into a product of unintegrated 
distributions, because the quark and the anti-quark can separately interact with 
the nucleus after they arc produced. The same is true also for the single quark 
distribution. In the limit of very large mass pairs, or when the quark and the 
anti-quark have large transverse momenta, fc_L-factorization is recovered. 

Though not fc^-factorizablc in the usual sense, the pair production cross- 
section can still be expressed in terms of the product of the unintegrated gluon 
distribution of the proton times the sum of four terms containing four such un- 
integrated distributions describing the nucleus. One of these four terms is the 
2-point correlator of two adjoint Wilson lines which already appeared in gluon 
production in pA collisions. This distributions corresponds to a gluon interact- 
ing with the nucleus in the amplitude and in the complex conjugate amplitude. 
Two of the other terms contain a 3-point correlator of two fundamental Wilson 
lines and an adjoint Wilson line. This corresponds to a gluon in the amplitude 
and a quark anti-quark pair in the complex conjugate amplitude or vice versa. 
The fourth unintegrated distribution is related to the 4-point correlator of four 
fundamental Wilson lines; this of course corresponds to quark-anti-quark pairs 
in both the amplitude and the complex conjugate amplitude. 

The single quark distribution is also not A;_L-factorizable. A significant sim- 
plification nevertheless occurs because the correlator of 4 fundamental Wilson 
lines is replaced by the correlator of two fundamental Wilson lines. One there- 
fore has only contributions from 2-point and 3-point correlators of Wilson lines 
to the single quark distribution. 

The three and four point correlators that appear in pair prodiiction cross- 
section probe non-trivial multi-parton correlations in addition to those involved 
in inclusive gluon production alone. The three point function that appears 
here is the same as one that appears in inelastic and diffractive deeply inelastic 
scattering at high energies [45,46]. These correlators of Wilson lines are therefore 
the relevant universal degrees of freedom at high energies, appearing in all high 
energy scattering reactions. The quantum evolution of these N-point correlators 
of Wilson lines at small x is described by the JIMWLK equations. Measurements 
of the pair production and single quark distributions can therefore provide a 
sensitive test of such evolution equations. 

In the McLerran-Venugopalan model of Gaussian correlations, these 2,3 and 
4-point correlators of Wilson lines can be computed analj^tically. The explicit 
expressions for these are derived in appendix B. The 2-point function has been 
known for a long time [31,47]. The 3-point function is, for Gaussian correlations, 
the g'^g-propagator and was computed in refs. [46,45,48,49]. For single quark 
production in the MV-model, as discussed here and by Tuchin previously, it is 
only these correlators that appear in the cross-section. For pair production, one 
needs the 4-point function as well. This is computed for Gaussian correlations 
in appendix B. 

As we discussed extensively in I, the MV-model is a good model for high- 
energy scattering when x is not too small - or quantum evolution is not im- 
portant. When quantum evolution sets in, the initial condition of Gaussian 
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correlations changes quickly [50,51] with non-Gaussian correlations becoming 
increasingly important at smaller values of x. The data at RHIC suggest that 
the MV-model might provide a good description of Deuteron-Gold collisions 
at central rapidities (when x ^ 10~^) but is a bad description already at for- 
ward rapidities corresponding to a; ~ 10""^. Computations of single quark and 
pair-production cross-section in both the MV-model and with the JIMWLK 
evolution equations will greatly help confirm the picture of D-Au collisions that 
is emerging from the RHIC data. We plan to pursue the phenomenological 
implications of our results in a future publication. 
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A Gaussian averages of Wilson lines 
A.l Introduction 

In this appendix, we derive the expression of the correlators of Wilson lines 
that appear in eq. (39). The most complicated one is the 4-point function 
C{x^,y^]u^,Vi_) = trc{U{xi_)t''W{yj^)U{ui_)t''W{vjL)), and it is trivial to 
verify that the 3-point and 2-point correlators that also appear in eq. (39) can 
be obtained as limits of the 4-point function: 



trc (uix±)t''U\y^)t^U^''\u±)) = C{x±, y^; u±,u±) 

tr^ (t''U'"'{xi_)t''' U'""''' {ui_)'^ =C(a?_L,a;±;M±,M±) ■ (55) 

Thanks to these relations, it is sufficient to compute the 4-point function. Simi- 
lar correlators have been studied in [46] and [52] , although the 4-point function 
wc need here is not derived explicitly in these papers. It is in general impossible 
to obtain closed expressions for these correlators. An exception is the case of a 
Gaussian weight function for the probability distribution of the source : 



WAPa] = exp 



J ^l^l{x+) 



(56) 



We denote-^^ P2{x~^,x±) = Se{x~^)p^{x±), and /Lt^(a;+) is the density of color 
charges at a given a;+ (this function is strongly peaked around a;+ =0). The 
most elementary correlator is 

(PA,a(a;+,a;_L)p^,6(y+,y_L)) = 6abS{x+ - y+)6{x^- yj_)ii\{x+) , (57) 



i^We must reintroduce the smearing in a;+ here in order to make sense of the path ordering 
that appears in the Wilson lines. It will however not play any role in the final results. 
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which is local in a;+ and x±. (The latter property is not essential in the following 
calculation, and can be relaxed if necessary.) 

In the calculation of the 4-point function C{x±,y^;u±,v±), the most com- 
plicated task is to deal with the SU{N) color algebra. A very convenient way 
to do that is to use a semi-graphical method. We can represent pictorially the 
correlator that we need as 



c 



->•• 



C{x_L,yj_;u_L,v±) = V— 4 . (58) 



c 



•>■• 




In this diagrammatic representation, the dotted lines represent Wilson lines in 
the fundamental representation with the horizontal axis the axis. Note that 
they therefore depend on the source to all orders. The arrows indicate the 
direction of the path-ordering, i.e. whether we have aU oiaU^ . The thick black 
dots are color matrices t° in the fundamental representation, and the wavy line 
that connects them denotes a contraction of their color indices (t" • • • i"). The 
Wilson lines are connected at a;"*" = ±00 by solid lines which merely indicate 
how the Wilson lines are multiplied (these solid lines do not contain any 
nor color matrices) . The closed loop indicates that the product of Wilson lines 
is traced. An essential ingredient in the forthcoming calculation is the Fierz 
identity in SUiN^^: 



ttA = - ^^iM . (59) 



which has the following graphical representation: 

Thus, the 4-point correlator C simplifies into a sum of two simpler terms: 




It is therefore sufficient to calculate the following object of SU{N) x SU{N): 

M^,kl{x±,y^,u^,v^) = (\u\v^)U{x^)] \u\y^)U{u±)\ \, (62) 

\ L Uj L ikl / 



^^See [53] for a detailed discussion of group algebra manipulations in non-abelian gauge 
theories. 
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which can be represented as a diagram with open ends: 



■Mijki{x±,yj_,uj_,v_i) = ^ . (63) 



Once we know this object, the correlator C is given by: 

'ikkiix±,y_i_,u±,v 

(64) 



C{x±,yj_;u±,v±) = ^Miikk{x±,y_i_,u±,v±) - -^Mikki{x±,yj_,u±,vj_) . 



A. 2 Tadpole corrections 

Thanks to the Gaussian form of the functional VK[p^], we can expand M. in 
powers of the elementary correlator {Pj^p^) using Wick's theorem. This amounts 
to expanding the Wilson lines in powers of the source and to connecting the 
Pj^ 's pairwise. From the point of view of its color structure, such a pairing {p^Pa) 
contains a contracted pair t'' ■ ■ - f^ of color matrices and can be represented by 
two black dots connected by a wavy link. Since {p^p^) is local in x+, only 
links that have their endpoints at the same x+ arc allowed. Some of them are 
"tadpoles" that start and end on the same Wilson line, and can be represented 
as follows: 



c 



(65) 



c 



In this diagrammatic representation, the black dots are color matrices^* 
and the solid line simply indicates in what order they are multiplied. A tadpole 
insertion on the Wilson line of transverse coordinate x_l at the "time" z'^ brings 
a factor —C pH'^^{z^)L{x,x), with: 

L{x,x) = g^ [ Go{x^-zi_f, (66) 



where Cp = (iV^ — 1)/2N is the quadratic Casimir in the fundamental repre- 
sentation of SU {N) and where 



(2^ feT 



Goixs. -z^)= I -j—^ (67) 



^®For tadpoles, there is one t" at each end of the wavy line, even if only one dot appears in 
the diagram. 
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is the 2-dimcnsional free massless propagator. Note that tadpole insertions 
are color singlet, and therefore commute with the rest of the expression. As 
a consequence, we can easily factor them out and resum them. This gives a 
factor^^: 

where we denote the density integrated over the longitudinal coordinate: 



+00 

,2 _ 



J dz+n\{z+), (69) 



by the same symbol, but without an argument. 

A. 3 Non-tadpole corrections 

We can now write 

M=TN , (70) 

where N is the same as M, but containing only links that connect different 
Wilson lines. A typical term in M with two links could be the following: 



c 



Note that since (p^Pa) local in 3;+, links must connect points at the same 
a;+. Therefore, diagrams with crossed links are not allowed. By a systematic 
use of the Fierz identity of eq. (60), we can completely remove the links. On 
the previous example, this would lead to the following four terms^°: 



(72) 



Note that the closed loop that appears in the second term can be cvahiatcd 
directly and replaced by the trace of the unit matrix in the fundamental repre- 
sentation trc(lj^) = N. Therefore, we see that all the terms that are generated 

^^The factor 1/2 in the exponential is due to the fact that the two ends of a Unk axe ordered 
in z+ if they belong to the same Wilson line. 

^"The diagrams have been rearranged a bit in order to make them look nicer, without 
changing their topology. 
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in the process of eliminating the Hnks with the Fierz identity are one of the 
following two terms: 



cz 



V 

a 



(73) 



This property, observed for the example of eq. (71), is in fact completely general. 

Let us now denote Afniz^, ■ ■ ■ , z+) the piece of Af that contains n links, at 
the "times" < ■ ■■ < z:^. We therefore have. 



J\f 



= E / M^(4,•••,4)• 

n=0 , / , 



(74) 



Using the above remarks, we can write each Mn as a linear combination of jV'^"^ 
and N^^'>: 

Afn = anM^"^ + hnAf^''^ . (75) 

Note that the term with zero links is exactly A/'^"\ so that we have ao = 1 

and 6o ~ 0- Once we know A/'n-i, the next term in the expansion, A/'n, can be 
obtained by adding one link to A/'^"-' and A/"^^^ in all the possible ways: 



J^n = /i^(^+)a„_i 



L{x, v) 



c 



L{x,u) 



c 



+ 



c 



c 



L{x, v) 




L[x, u) 




, (76) 



where, generalizing eq. (66), we denote 

L{x, v)= g'^ Go{x_i_ - z±)Gq{v± - z±) . 



(77) 



There are six possible ways to link four lines pairwise, and we have only repre- 
sented two of them. Note that there is a — sign for links connecting two Wilson 
lines going in the same direction (i.e. two J7's or two JJ^'s). The next step is to 
use eq. (60) in order to remove all the links and to replace by the closed loops 
that may appear in the process. Finally, identifying with A/'„ — a„A^("' +6„A/'*^''' , 
we obtain the following recursion relation: 



(78) 
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where M is a 2 x 2 matrix given by 



(79) 



with 



a = L{x,v) + L{u,y) , 
(} = L{x,u) + L{v,y) , 
^ = L{x,y) + L{v,u) . 



(80) 



Note that M does not depend on z+. Given the initial condition ao = 1, 60 = 0, 
it is now very easy to formally solve the recursion relation to obtain 



i=l 



(81) 



From the coefficients an,bn, the correlator C{xj_,y^;uj_,vj_) is then ob- 
tained by using eqs. (61), (70), (74) and (75). Using eq. (61) amounts to "closing 
the ends" of Af. Multiplying also by the contribution of the tadpole insertions, 
we get: 

+ OC 



C{x±,yj_;u±,v±) = ^ 



E 



(iV2 - 1) 



T 



(82) 



In order to make this formula more explicit, we must integrate eq. (81) over the 
z^'^s, which gives^^: 



1 



n! 



(83) 



z:l;<---<zt 



The only remaining task is to calculate the n-th power of the constant matrix 
M. This can be done by finding its eigenvalues, which are easily shown to be 



A± = 



1 



N 1 

Y(a + 7) + ;^(/3-a-7)± 



iV2 



(«-7)' + (a-/?)(7-/3) 



It is then clear that we can write 

+00 



E 

n=0 



(84) 



(85) 



z+<---<zt 



^^The ordering of the z^'s in eq. (74) can be trivially removed because the recursion matrix 
commutes with itself at different times. This gives a factor 1/n!. 
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with a-|_ and a_ two numbers determined from the initial condition. At this 
point, it is a simple matter of algebra to derive our final expression for the 
4-point function C: 



xe-(^-4Tv)M^ [r{x-v)+r(x-y)+r(u-v)+r(u^y)] 
I 2\/A 2\/A I 



where we denote 



and 



1 



16 (a-/3)(7-/3) 



{a - 7)2 

r(y -u) = L{y, y) + L{u, u) - 2L{y, u) . 



(86) 
(87) 
(88) 

A comment is in order regarding this result. The quantity L{x,y) defined in 
eq. (77) is strongly infrared divergent. Indeed, if one introduces by hand an 
upper cutoff in the integration over z± at z± ~ J^qco^ then L{x,y) ~ ^qcd- 
However, the combination T of these objects, defined in cq. (88), has a much 
milder infrared singularity which is logarithmic in the cutoff scale. This is also 
the case of the differences a — 7, a — /3 and 7-/3 which appear in eq. (86). 
Therefore, we conclude that the 4-point function studied here is not more in- 
frared singular than the 2-point functions already discussed extensively in the 
literature. Physically, the reason why the quadratic divergence is softened to a 
logarithmic divergence is that we consider a correlator made of a singlet com- 
bination of Wilson lines. The reason why we have nevertheless a residual log- 
arithmic divergence is that the color neutrality of the nucleus is not enforced 
in the functional of eq. (56) . We have introduced by hand a cutoff at the scale 
^QCD which is equivalent to imposing color neutralization at the scale of the 
size of the nucleon. 



A. 4 Large limit 

Eq. (86) becomes substantially simpler in the large N limit. Indeed, in this 
limit, we have A = 1. We have therefore, 

iV— >(x) 2 

(89) 

Using the explicit expressions of a, 7 (eq. (80)), this simplifies greatly to read, 
C{x^,yr,n^,v^) = ^e-*''lP(^-^)+^(^-")l. (90) 

JV— >oo Z 

This result could also have been obtained very simply from the recursion relation 
eq. (78) because the recursion matrix M is diagonal in the large N limit. 
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A. 5 Special cases 

There arc a few degenerate cases in which the 4-point function we have calcu- 
lated in this appendix simplifies into a simpler function. These special cases are 
useful as a check, because for these the result is already known in the literature. 

A. 5.1 x± = v± or = u± 

There is an obvious simplification in C{x±,yj_;u±,v±) if x± = v± or if j/j^ = 
u±. Indeed, any of these equalities makes /3 = 7, which implies A = 1. For 
instance, in the case of x± = v±, we obtain from eq. (86) the following formula: 

Cix^,y^;u^,x^) = iVC^e-^^^A 

= TVC^e-^^^^"-"^ . (91) 

This simple formula is indeed what we already know (see [39] for instance) for 
this degenerate case, since we expect: 

C{x^,y^;u^,x^) = C,ti:,(u\y^)U{u^)) . (92) 

A. 5. 2 x± = y^ or u±^ = v± 

A slightly less trivial simplification is obtained when x± = y^ or uj^ = v±. 
For instance if u± = vj_, the correlator C should become the 3-point function 
trc{U {x±)t°'U^y 2_)t''U''"- {u±)'^ , whose expression can be found in the appendix 
of a paper by Kovner and Wiedemann [46] . In the c;asc where one of these two 
equalities is satisfied, we have a = (3, which also implies A = 1. When u± = v±, 
eq. (86) simplifies into: 

C{x^,yr,u^,u^) = iVC^e-wA'l [r(-«)+r(y-«)] 

Xe-(f-2T^)^l [r(x-«)+r(a;-j;)+r(«-y)]gf M^(a-7) 

= NC^ e~^^'^ [rix-u)+r{y-u)]-[^-f)f.lr{x-y) ^^^^ 

It is trivial to check that this formula is equivalent to eq. (C7) of [46] (the 
correspondence between our T{x—y) and the v{x—y) defined in [46] is 2v{x—y) = 

B Retarded pair production amplitude 

B. l Introduction 

In section 3, we have calculated the time-ordered pair production amplitude, 
which is directly related to the probability Pi of producing a single quark- 
antiquark pair in a collision. We also conjectured that, since we are treating 



29 



the problem at leading order in the source pp that describes the proton, at most 
one pair can be produced in a collision. This implies that the average number 
of pairs Nqq produced in a collision must be equal to the probability Pi , simply 
because the probabilities P2, P^,- " to produce 2 or more pairs in a collision 
are all zero at this order. 

On the other hand, we can calculate directly the average pair multiplicity 
Nqq from the retarded pair production amplitude, as explained in [38] , from the 
formula, 

where M.j^{q,p) is the pair production amplitude calculated with retarded 
quark propagators (as opposed to Feynman quark propagators in the case of 
Ai p{q;P))- Since the calculation of the retarded amplitude presents a number 
of significant differences with respect to the time-ordered amplitude, its calcula- 
tion provides a non-trivial cross-check of the results we have derived in section 
3. Because of the expected equality Pi = Nqq, we expect the retarded and 
time-ordered amplitude to be identical "up to a phase" . 



time-ordered : 



retarded : 




Figure 3: Top: the contributions to the Feynman amplitude. The shaded area 
represents the nucleus. A gluon is emitted by the proton, and splits into a 
pair: after the collisions with the nucleus (left), inside the nucleus (middle) or 
before the collision (right). Bottom: the same contributions for the retarded 
amplitude. 

In fact, without doing any calculation, we can guess what this phase is. The 
contributions to the two amplitudes are sketched in figure 3, and one can see 
that the only difference between the two is that the antiquark line comes from 
the region of negative in the case of the retarded amplitude. This is so 
because the retarded quark propagator can only move forward in time. From 
this figure, we see that moving the antiquark line from the left of the diagram 
to the right of the diagram changes the amplitude by a factor which is precisely 
the scattering matrix of the antiquark on the nucleus. In a sense to be made 
more precise later, we expect the following relation: 

= M, U{y^) , (95) 
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where yj^ is the transverse coordinate of the antiquark. 

In order to calculate explicitly the retarded amplitude, we must use the free 
retarded quark propagator: 

SliP)^i , ^tT- + ^ (96) 
" — m"' + ip+e 

which differs from eq. (18) by the location of its poles in the complex energy 
plane. Note that the p'^e in the denominator could as well have been p~e or 
p^e because p'^,p~ and have the same sign at the pole. The other technical 
difference is the form of the scattering matrix of a quark/ antiquark on the 
nucleus. In the retarded case, it is given by [38]: 




27r5(fc+)7+ J (fx±e'''^-''^ [u{x±) - l] . (97) 



We see that it coincides with the time-ordered scattering operator given in 
eq. (16) when p+ > (as expected since the poles of the retarded and time- 
ordered free propagator are the same if p+ > 0), but differs from it whenp+ < 0. 

B.2 Singular terms 

Let us start with the calculation of the term involving the field ^^i„g, where the 
pair is produced inside the nucleus. It is trivial to reproduce the calculation 
of section 3.3 with the Feynman rules appropriate for the retarded propagator. 
We obtain: 



This is indeed the same as eq. (31), up to the factor U{x±) (the transverse 
coordinate of the antiquark is x± in this expression). 



B.3 Regular terms 

We must now calculate the four "regular" (in the sense that the pair is produced 
outside of the nucleus) diagrams of figure 1 in the case of the retarded amplitude. 
The diagram (a) does not contain any internal quark propagator and is therefore 
identical in the time-ordered and retarded cases: 

M^^\q,p)=Mi^\q,p) . (99) 



31 



Contrary to the tinic-ordcrcd case, the pole structure of the diagram (d) is such 
that it does not contribute in the retarded case^^: 

M'^;}\q,p) =Q . (100) 

The diagram (c) is also very simple: indeed, the only fermion propagator that 
interacts with the nucleus is the quark line which carries the momentum g+ > 0. 
Therefore, this propagator is the same in the retarded and in the time-ordered 

case: 

M^^\q,p) = M(^\q,p) . (101) 

All the difficulties are concentrated in the diagram (b). Indeed, since this di- 
agram has a fermion propagator carrying the momentum — p+ < 0, its pole 
structure is not the same in the retarded and in the time-ordered case. In addi- 
tion, the interaction operator of the antiquark on the nucleus is different. The 
diagram (b) can be divided in three terms: 

= + M% + ^«,3 • (102) 

■M^^\ is obtained by keeping only the first line of the field in eq. (5) (namely 
replacing the solution of Yang-Mills equations by the field of the proton alone). 
This term can be written as follows: 

u{q)l+{i-t + m)-i-{i-$-ti + m)7+[i'^-t''C/(y^)];;(p) 



2p+((g± - + to2) + 2q+((q^ - k± - ki±)'^ + m"^) 



(103) 



The term we denote ■M.'"^\ is obtained from the diagram (b) by keeping only 
the piece proportional to C^ ,.^g in eq. (5). Using the fact that 

Ci^,re,{P + q-k) = 2{p + q-kr- 5^- ^l^^ X+ ' ^^""^^ 

we obtain: 

dk- u{q)t'>0^,rcs{p + q - k)i-l^ + ^ + m)-f+[Uix±) - l]v{p) 

2w [{p + q- fc)2 +i(p+ + q+ - k+)e] [{-p + kf - + i{-p+ + k+)e] 

^ . uiq^+t'lUix^) - Ijvjp) 

p+ + q+ ■ ^ ' 

With the help of this relation, we obtain easily: 

.Jb) _ 2 f Pp,a(fci_L) 

x?«2;^%4^[n«x)-ii.. (106) 



f d ki^pp^ajki^) r ,2 ,(p^+q^-fei^).a,^ 

J {2^r kl^ J 



^^This property is a consequence of causality. Indeed, in the retarded case, the diagram (d) 
corresponds to the following process: a quark goes through the nucleus, absorbs the gluon 
coming from the proton, and goes again through the nucleus. At high collision energy, this is 
suppressed by one power of the collision energy. 
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Finally, -M^^"^^ is obtained from the diagram (b) by keeping only the piece pro- 
portional to in cq. (5). It reads: 

~ 7 (27r)2 (2^)2 kj^ jax^ay^e e 
dk- u{q)<p^{p+q-kMA){t-I^+rn)-(+t\U{v^)-l\v{p) . 



27r [ip + q-kf+ieWi-p+kf-m^-ie] "[^('^^^ 

(107) 



B.4 Total retarded amplitude 

The main difficulty at this point is that wc must combine this term in (p + 
q — k, ki±) with a term in C^{p + q, ki±) contained in In order to do 

this, one must notice the following relation: 



u{q)t^^„{p + q,ki±)v{p) 



'I 



{p + qY 

dk- u{q)t''(p„ (p + q-k, fci_L)(-2/ +^ + m)7+w(p) 



277 [{p + q-kf + ie][{-p+kf 



(,108) 



which is valid if fex = in the r.h.s. and introduce in A^^"'' dummy variables 
y_i_ and fe_L via a trivial factor (which ensures that A;_l = 0): 

1 = / '^.d'v^-''-''- ■ (109) 

In order to simplify the notations, let us denote: 

dk- (:^,{p + q~k,kn_){-j/+l^ + m) 



^{q±,P± - kj_,ki±) = i 



27r [{p + q- + ie] [{-p + k^ - - ie] ' 

(110) 

where we have only indicated the transverse momenta in the list of arguments. 
Note that this object depends only on the difference pj_ — k±_ rather than pj_ and 
A;_L separately. Finally, we need to perform the transformation fe j_ — » Px + Q± ~ 
fcix — fc± for the terms containing C^j (or equivalently the T defined above). 
Combining all the contributions to the retarded pair production amplitude, we 
can write it as follows: 

MAq,p) = J d^y^e'P^ y^u{q)IC{qr,yjU{y^h+v{p) , (111) 
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with 



A-2 



^2jj.^gife± -a: J. ^i{q^ -k± -ki±)-yj_ 



(27r)2 (27r)2 



+ 



t''^{Qj_, k± + ki± -q±,ki±)- 



p+ +q+ 



[U{x±) - l]ba } . (112) 



Note that the quantity )C{q±; y±) does not depend on pj^. 

B.5 Compcirison with the time-ordered amplitude 

Let us now go back to our original problem, namely to verify that we indeed have 
Pi = Nqq. In order to make this verification easy, the timc-ordcred amplitude 
given in eq. (34) needs some rewriting. By performing on the term in C„ the 
transformations described in the section B.4, we can write: 

MAq,P) = J d'y^e'^^-''^u{q)JCiqr,y^h+v{p) , (113) 

with the same Kl{qj_; yj_) as in eq. (112). We see that the time-ordered and the 

retarded amplitudes differ only by the unitary matrix U{y^) in the integrand: 
this is the precise meaning of "identical up to a phase" in this context. 
From eq. (113), we can write the single quark spectrum as: 



dPi 



1 f (Pp^ dp+ f ^ 2 / 
4^^ J J2^^ J ' 



^p±-(y±-yi 



d'^q^dyq 6477* 

xtrd[(^ + m)/C(q^;y^)7+(2/-m)7V/Ct(q^,y'^)7°] • (114) 

Since m)j~^ = 2p+7+ does not depend on nor docs /C, it is trivial to 

perform the integration over pj^, which gives a S{yj_ — y'j_). Therefore, we can 
rewrite the above equation as follows: 



dPi 



i-j j dp+d^y^tT^[{4+m)K:{qr,yi.h+l°}C\q^,y^h°] . (115) 



d'^q^dyq 3277* 

Similarly, one can write the average number of produced quarks as: 
dN„ 



I qq 



1 / d^p^ dp+ [ 2 ,2 / 



ip±-iy±-y'±) 



d^q^dyq 64n^ 

xtra [{4 + m)lC{q^; y^)U{y^h+(^- mh+W {y'^)^°IC'' [q^Mj^] . 

(116) 

Again, the integration over pj^ produces a 6{yj^ — y'±)j and we see that the ma- 
trices U{yj_) and U^{y'j_) cancel each other to produce a unit matrix. Therefore, 
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we conclude that: 

rfPl _ dNgg ^^^^^ 



<Pq±_dyq (Pqj_dyq 
as expected. 
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